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Abstract: Riemannian cubics are curves in a manifold M that satisfy a variational condition appropriate for 
interpolation problems. When M is the rotation group SO(3), Riemannian cubics are track-summands of Rieman- 
nian cubic splines, used for motion planning of rigid bodies. Partial integrability results are known for Riemannian 
cubics, and the asymptotics of Riemannian cubics in 50(3) are reasonably well understood. The mathematical 
properties and medium-term behaviour of Riemannian cubics in SO (3) are known to be be extremely rich, but there 
are numerical methods for calculating Riemannian cubic splines in practice. What is missing is an understanding of 
the short-term behaviour of Riemannian cubics, and it is this that is important for applications. The present paper 
fills this gap by deriving approximations to nearly geodesic Riemannian cubics in terms of elementary functions. 
The high quality of these approximations depends on mathematical results that are specific to Riemannian cubics. 
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1 Introduction 

Suppose that a C°° curve x : K — > M in a Riemannian manifold M is sampled at times to < t\ < ... < t n , 
yielding observations x(ti) — Xi for < i < n. Then x is uniformly approximated by a track-sum of minimal 
geodesic arcs joining successive observations, with 0(5 2 ) error, where 5 is the maximum distance between Xi-\ and 
Xi, for 1 < i < n. So any C°° curve in M is a track-sum of curves that are nearly geodesic. Although x is C°°, the 
piecewise-geodesic approximation usually fails to be C 1 at junctions. 

A C 2 approximation is given by a natural Riemannian cubic spline, namely a track-sum of Riemannian cubics, 
critical for the mean squared norm of covariant acceleration. If sampling is sufficiently frequent we can restrict 
attention to Riemannian cubics that are nearly geodesic. 

In the special case where M is flat, Riemannian cubics are expressed simply in terms of cubic polynomials. When 
M is curved, Riemannian cubics are given by a differential equation that is difficult to solve, even when M is the 
unit 3-sphere with the standard metric, or the bi-invariant rotation group 50(3). These cases occur in motion 
planning for rigid bodies, and so numerical methods are needed to find Riemannian cubic splines 0. 

The asymptotics of Riemannian cubics are studied in [TO], [TT], [T3], but little is known about the short-term 
behaviour of cubics, and the short term behaviour is more relevant for motion planning of rigid bodies. The present 
paper fills this gap, by deriving approximations to nearly geodesic Riemannan cubics in terms of trigonometric 
functions and polynomials. The new approximations are much more informative than Taylor approximation, either 
in coordinate charts or ambient space, and capture interesting short term geometry of Riemannian cubics that was 
previously observed in numerical experiments. 

Before describing the methods and layout of the present paper, we review Riemannian cubics in more detail. 

2 Riemannian cubics 

For M a finite-dimensional Riemannian manifold, consider the functional 

J(x) := f\v t x^(t),V t x^(t)) dt 
J t 

defined on C°° curves x : [to,ii] — > M, where x and its derivative x^ are prescribed at to,t\. Here V denotes the 
Levi-Civita covariant derivative defined by the Riemannian metric ( , ). A Riemannian cubic is a critical point x 
of J, in the sense that x satisfies the associated 4th order Euler-Lagrange equation [9] 

V?x«+i^ (t) (V t z« = (1) 

where R is the Riemannian curvature. The existence of Riemannian cubics satisfying the prescribed conditons is 
proved in [5] when M is a complete Riemannian manifold. As seen from ([!]), cubically reparameterised geodesies 
are Riemannian cubics. However most Riemannian cubics do not arise in this way. 

Definition 1 Given e > 0, a C°° curve x : \to,t\\ — > M is (e-Jnearly geodesic when, for all t £ [to,ti], 

||V t a; (1) (i)|| < e and ||V t 2 a; (1) (<)|| < e. 

□ 

Riemannian cubics need not be nearly geodesic, but the restriction of any C°° curve to a sufficiently small subin- 
terval can be reparameterized to a nearly geodesic curve defined over a fixed interval [£o, £i] ■ So nearly geodesic 
curves are informative about the local geometry of arbitrary C°° curves, in particular Riemannian cubics. Nearly 
geodesic curves arise naturally in other ways too. 

Example 1 Take M to be the matrix group SO(3) of rotations of Euclidean 3-space E 3 . Define x(t) £ 5*0(3) 
by taking its columns to be the coordinates at time t of an orthonormal frame fixed relative to some rigid body B. 
If the mass distribution of B is spherically symmetric, and if B moves freely, then x is a geodesic in 50(3) with 
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respect to a bi-invariant Riemannian metric. If, however, B is subject to a C 1 -uniformly small torque T then x is 
nearly geodesic. 

Let T be unknown, and suppose that the configuration of B and its angular velocity are observed at times t ,ti. 
Then a minimis er x of J is an interpolant that minimises the mean-squared torque. Since T is C l -uniformly small, 
the interpolant x is a nearly geodesic Riemannian cubic. □ 

We refer to [H [TJ1 [TSl GH H] for further applications of Riemannian cubics. 

When M is Euclidean m-space E m , a Riemannian cubic is precisely a polynomial curve of degree < 3. Nothing like 
this can be said for non-flat manifolds M, even when M is a space of constant nonzero curvature. The situation 
for elastic curves is entirely different. 

Example 2 An elastic curve in M is a critical point of the restriction of J to the space of constant-speed curves 
x, with x and x^- 1 ' still prescribed at to,ti. When M — E 3 elastic curves are the Euler elastica, whose curvature 
and torsion are obtained in terms of elliptic sine function, as in Lecture 1 of [18j . whereas Riemannian cubics in 
E 3 are just cubic polynomial curves. 

On the other hand, when M is the unit sphere S 3 in E , there are quadrature formulae for elastic curves in terms 
of elliptic functions [H [1] I17j and Lecture 2 of |18) , but quadrature formulae for Riemannian cubics in S 3 are known 
only for a codimension 3 subclass [16] . Elastic curves do not resemble Riemannian cubics, except for the very short 
term. □ 

Whereas the long term behaviour of Riemannian cubics is studied in |10j , |12j , , their short and medium term 
behaviour is poorly understood. Yet the short and medium term are more significant in applications, such as 
interpolation and motion planning for rigid bodies. 

Example 3 Given to < t\ < . . . < t n and Xq, X\, . . . , X n £ M define 

J(x) := f n (V t x<»(t),V t xV>(t)) dt 

defined on C 2 curves x : [to, t n ] — > M satisfying x(ti) = x. L for % = 0, 1, . . . ,n. A critical point of J is called 
a natural cubic spline. Natural cubic splines are characterised as C 2 track sums of Riemannian cubics on the 
intervals [tj_ i, ti] for i — 1,2, ... ,n, whose covariant acceleration vanishes at to and t n . □ 

As we shall see in Example [9j Taylor approximations are of limited value in this context. Much more accurate 
estimates can be made by exploiting specific properties of Riemannian cubics x, especially when x is nearly geodesic. 

We focus on bi-invariant M = SO(3) of rotations of E 3 , and on Lie quadratics namely the left Lie reductions 
V : [^Oi^i] — ^ E 3 of a Riemannian cubics x. The Lie quadratic Vs of a nearly geodesic Riemannian cubic xg is 
nearly constant. The variational equations of Lie quadratics are used to find approximations of Vs. There is a 
quadrature formulae of [12] for Riemannian cubics in terms of Lie quadratics, but the approximate Lie quadratics 
cannot be directly substituted for Vs into the formula. Nonetheless, taking a little more care, we obtain an 
approximation x to xs- 

In effect, the known structural results for Riemannian cubics in S*0(3) are exhausted, before resorting to a Taylor 
approximation in what is left. Surprisingly whereas the reconstruction of Riemannian cubics from Lie quadratics 
requires a quadrature [12j . the first order approximation x for xs given in Theorem [4] is algebraic]^] in trigonometric 
functions and low degree polynomials. 

The layout of the paper is as follows. 

• [|3] is an introduction to Lie reductions of Riemannian cubics in bi-invariant Lie groups G. 

• Q studies the variational equation of a Lie quadratic, giving examples where the variational equation can be 
solved exactly. Example [6] concerns the simplest case of variations through nearly constant Lie quadratics. 
Derivatives V"^' ) to order < i < n of a variation to order n with respect to the variation parameter, give 
rise to an order n approximation V n of a nearly constant Lie quadratic Vs. 

The expression is complicated, but it is difficult to see any way around this. 
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• In Sj5j C? is taken to be the Lie algebra so(3) = E 3 of the rotation group 50(3). Theorems [2J [3] give formulae 
for the first and second order approximations V\ and V2 to the Lie quadratic V§. In applications such as 
Example [T] Vg gives the angular momentum relative to the body B. 

Already V\ captures significant medium term behaviour of Vg, as illustrated in Figures [lj [2] of Example [9} 
Taylor approximations to Vg perform badly (Figure [l]), while V\ and V% are nearly indistinguishable from Vs 
in the medium term. For the longer term, V% significantly improves on V\ (Figure EJ). 

• In sj(3] the approximations V n for n = 1 and n — 2 are used in combination with Theorem 5 of |12) to derive a 
first order approximation x to the Riemannian cubic xg, namely Xg(t) — x^\t) + 0(S 2 ) for any non- negative 
integer j. This is our main result, stated as Theorem [4j and implemented in Example |10| 



3 Riemannian cubics in Bi-Invariant Lie groups: Lie quadratics 



Now we take M to be a path-connected finite-dimensional Lie group G, with bi-invariant pseudo-Riemannian 
metric. The restriction of the metric to the Lie algebra Q is an ad-invariant semi-definite inner product, namely 
ad(it) : Q — > Q is skew-adjoint for all u G Q. Conversely an ad-invariant semi-definite inner product]^] on Q 
extends by left multiplication to a bi-invariant pseudo-Riemannian metric on G. Define the left Lie reduction V of 
x : [toM] -> G by 

V{t) :=dL(x(t)-% it) (x^(t) 
where L(g) denotes left multiplication by g € G. 

Theorem 1 ([9], llOp x is a Riemannian cubic in G if and only if, for all t G [to,ti], and some C G Q , we have 

V^(t) = [V^(t),V(t)] + C, (2) 

□ 

Equation ^ is second order, whereas the Euler-Lagrange equation ([T]) has order 4. Equivalently we may write 

1*3) (t) = [V^(t),V(t)}. (3) 

A curve V : [toi^i] ~~ ^ G satisfying |2]) for all t is said to be a Lie quadratic. The Lie quadratic and Riemannian 
cubic are said to be null when C = 0, and non-null otherwise. Null Lie quadratics in E 3 appear in applications in 
fluid dynamics [7]. 

If x : [to, t{\ — > G is a Riemannian cubic then so are x~ x and gx where g G G. Geodesies in bi-invariant Lie groups 
are precisely the Riemannian cubics with constant Lie quadratics. 

Example 4 An affine line in Q is a Lie quadratic, and is null when its image contains 0. A Lie algebra Q is 
abelian precisely when all its affine lines are null Lie quadratics. Then all nonconstant null Lie quadratics are 
affine lines, and a Lie quadratic is precisely a polynomial curve in Q of degree at most 2. □ 

From Theorem Q] a short calculation shows 

Corollary 1 Let y,z : [to, ^l] G be restrictions of 1-parameter subgroups of G, with infinitesimal generators 
A, B G Q . Then the pointwise product 1 1— > y(t)z(t) is a Riemannian cubic if and only if [[[B, A], A], B] =0. □ 

Example 5 Let G be bi-invariant SO(3), with A, B G so(3) linearly independent. Then the pointwise product 
1 1 y cxp(tA)cxTp(tB) is a Riemannian cubic if and only if (A, B) = 0. □ 

Null Lie quadratics in E 3 are studied in [101 114) . The present paper focuses on the non-null Lie quadratics and 
Riemannian cubics, which are generic, and whose behaviour can be complicated, as illustrated in |llj . 

2 These do not exist for some Lie groups. However for G semisimple we may use the Killing form. 
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4 Variational Equations For Lie Quadratics 



For t < t 1 and any C°° map V : [—1,1] x [t , *i] -> 5, define T4 : [£ , ii] ->• £ by 14 (t) = V(h,t), where 
ft e [—1, 1]. Denote the i-fold derivative of V with respect to ft of the j-fold derivative with respect to t, by V^'^. 

Suppose that, for each ft £ [— 1, 1], 14 is a Lie quadratic. Differentiating V^ ' 3 - 1 = [y(°< 2 ) ; l/(°> )] n > 1 times with 
respect to ft., we obtain the third order linear ODE for Y„(t) := V^ 11 ' ' (h,t) 

+ ad(^°))yi 2 ) -ad(y(°- 2 ))y n - E ( " ) ^,^'1 (4) 

i=l ^ ' 

where the coefficients and right hand side are in terms of derivatives with respect to t of Yi with i = 0, 1, 2, . . . , n— 1. 

Let Vb be the Lie quadratic of a nontrivial cubically reparameterised geodesic, namely Vq = qo(t)D for ^ D E Q 
and qo : K — > K is a quadratic polynomial. For n = 1, the variational equation Q evaluated at ft = becomes 

YF\t) + a A(D)(q (t)Y{ 2 \t)-q< i ) Y 1 (t)) = 0. (5) 

Let ad(-D) : Q — > Q be diagonalizabl^] over C, namely Q = K, + £ where K, is the kernel of ad(D), and £ has a basis 
of eigenvectors E^ of ad(D). For K € K, (|5j| gives 

(if^),^) = -(ad(D)(q (t)Y£ 2 \t) - <?< 2 Vi(t)), K) = (q (t)Y} 2 \t) - q {2) *i(f), ad(D)tf) = 0. 

So (Y±(t),K) = quit) where qK '■ \to,t\\ — > R is another quadratic polynomial. Similarly, (Yi(t),Ei) satisfies 

y^(t) = r(t)y^(t)-r^y(t) =► j/ 2 )(t) - r(%W(i) - r«(%(t) + Co (6) 

where r(£) := Ajgo(i), A, £ C is the eigenvalue of of Ei, and cq is constant. The linear ODE ([6| is solvable by 
quadratures in terms of solutions of the associated homogeneous ODE 

y^{t) = r(t)yW(t)-rW(t)y(t). (7) 



Example 6 Let r be a constant ^= k g C. TTien reads y^ 2 \t) — ky^\t) => y{t) = C2 + cie fct where c\,&2 
are constant. So the general solution of is 

y(t) = -|t + c 2 + c ie fe *. 

□ 

Example 7 Lei r &e nonconstant and linear. After a time-shift, write r(t) — kt where 7^ k G C is constant. 
Then /[ly becomes y^ 2 \t) = kty^(t) — fcy(t) and, substituting y(t) = tz(t), 

2z^+tz^ = kt 2 tz^ =>. = C2 ;-V* 2 / 2 

/rom which z(t) is found by quadrature. □ 

Example 8 Lei r = fci 2 where k is constant. Equation becomes ?/ 2 '(t) — kt 2 y^(t) — 2kty(t). Substituting 
y(t) — w(s) with s — kt 3 /3, gives a /orm 0/ Rummer's equation 

2 2 

sw = (s )w H — to 

v 3 ; 3 

where w' and w" are derivatives with respect to s. So w is given in terms of confluent hypergeometric functions. 
□ 



3 In a scmisimple Lie algebra |19| . the so-called semisimple elements D for which this holds comprise an open dense subset of Q. 
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Suppose now that Vb is constant, and let h i— > V(h,0) be an affine function of uniform norm < 1. From equation 
([6]), for any non-negative integers n and j, for some some c nj > 0, and for all (h,t) G [—1, 1] x [to>ii]) 

\\V^(h,t)\\ < Cnj. 

For 6 € (0,1) and n > 1, the ni/i order approximate quadratic V n : [to,ti] — > Q of V : [—1,1] x [t , ^i] — > Q is 
defined by 

i=l 

By Taylor's Theorem, Vj)?\t) = (f) + 0(5 n+1 ) where the asymptotic constants depend on Vq, n and j. Once 
Q is specified the analysis of Example [6] can be taken further. 

5 Approximating Nearly Constant Lie Quadratics in so(3) 

Considering Euclidean 3-space E 3 as a Lie algebra, with respect to the cross-product as Lie- bracket, the Eu- 
clidean inner product ( , ) is ad-invariant. The Lie isomorphism ad : E 3 — > so(3), given by ad(v)(u>) :— v x w, 
identifies E 3 with so(3). In particular, an ad-invariant inner product ( , ) on so(3) is defined by requiring ad to 
be an isometry. Take G — 50(3) with the corresponding bi-invariant Riemannian metric. 

We seek nth order approximate quadratics V„ : [to,ii] — > so(3) of variations V, where Vq is a nonzero constant 
D e so(3) and h i— > V(/i, 0) is affine. Then ad(Z?) is diagonalizable, the linear span H of D is a Cartan subalgebra, 
and ad(Z?) : so(3) ®r C — > so(3) ®r C has eigenvalues ±di, with unit eigenvectors E\,E% = E\ e so(3) <&g C 
orthogonal to D. Define unit vectors 

F :=D/d, Fx := (E 1 +E 2 )/V2, F 2 := - {E x - E 2 )/(V2i) 

with [Fo,Fi] = F2 and [i^-Fb] = ^l- Then Fo,Fi,F 2 corresponds under ad to a positively oriented orthonormal 
basis of E 3 . 

For n > 1 define /„ : [to, ti] — > K and u„ : [i , £1] — > Fq 1 - by f n (t)F + v n (t) := Y n (t). Taking k = ±di in Example 
|6j /1 is a quadratic polynomial q, and v\ — A + e(B) where A: [to, t{\ — > i 7 ^ C so(3) is affine, -B G so(3), and e(t) 
is the Lie endomorphism exp(— d(t — to)ad(Fo)) of so(3). We have proved 

Theorem 2 TTie /irsi order approximate quadratic of V has the form 

V x {i) = D + 5(q(t)F + Ao + (t-t Q )A 1 +e(t)(B)) 
where q(t) = cq + c\ (t — to) + c 2 (t — to) 2 , cq, c%, c 2 , d £ R and Aq, A\, B £ F^ . Then, for j > artrf t € [to, t\], 

V s u) (t) = vf j) {t) + 0(5 2 ) 
where the asymptotic constants depend only on d,to,t\ and j. □ 

For n > 2, integration of equation Q gives 

= £ ( J ) (8) 

^ - £ ( * ) i~eI(fi%e-\ - ffe-^) (9) 
i=i ^ ' 

where i := ad(i*b) and, for a continuous curve g of linear endomorpisms of so(3), I(g)(t) :— J tg g(s) ds. 
Then i commutes with e, and i 2 = — 1 where 1 is the identity on so(3). If g is C 1 , then 

/(e fl ) = ^(eg- ff (i )-/(g<7 (1) )), in particular 7(e) = ~(e-l). (10) 
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L := 




Li := 




M := 




Mi := 




M B := 





Taking n = 2 in equation Q, and writing j4(t) = An + (i — in)Ai with An, Ai € so(3), repeated use of (10) gives 
/a = -2<[>4o,io(B)] + L4i, L x (S)],F ) 

W 2 = 2q^(M Q (A ) + M 1 (A 1 )-M B (B))+2d 2 ioI 2 (I(q)e)(B) 
where Lo(t), L\{£), Mo(t), M\(t), Mb it) are the endomorphisms of so(3) given by 

(u 2 /2-l)i+<oe) (11) 

-3)l + 2ui + 3e + uiog) (12) 

-l)l + ui + e) (13) 

-u)l + {u 2 /2-l)i + ioe) (14) 

L) + uio(g+l)) (15) 

with e evaluated at t, and u :— d(t — to). This proves 
Theorem 3 The second order approximate quadratic of V is 

V 2 (t) = D + 5(q{t)Fo + Ao + (t-to)A 1 +e(t)(B)) + S ^{f 2 (t)Fo + v 2 (t)) 
where q, An, A±, B are as before. For j > and all t 6 [to>£i]; 

vp\t) = V 2 u \t) + 0(S 3 ) 

where the asymptotic constants depend only on d,to,t\ and j. □ 

Notice that q, A, B contain 3 + 4 + 2 scalar variables, sufficient for initial conditions of the 3rd order ODE ^ in 
so(3). 

Example 9 Using ad to identify so (3) with E 3 , Figure^plots (blue) a numerical solution of ^ for a Lie quadratic 
Vg : [0, 5] — > so(3) = E 3 near Vq = (1, 0, 0). The numerical solution is obtained using Mathematical NDSolve, and 
Vs is the non-null Lie quadratic specified by the initial conditions 

V s (0) = (1.005,0.006,-0.01), V^O) = (-0.005,-0.00449,0), 

V^ (2) (0) = (0.001,-0.005,0.005), V s (3) {0) = (0.00002,0.005035,0.005031) 
with C w (0.0009551,-0.00495,0.00051755). 

The approximations (green) V\ and (red) V 2 are nearly indistinguishable from Vs. The second degre^ Taylor 
polynomial, also shown (dashed) in Figure [IJ is a poor approximation to Vg. The initial dot labels Vs(0), and the 
second dots correspond to t = 2. 

Figure ^plots the same Vs (blue), V\ (green) andV 2 (red) fort 6 [0,25]. At first, both approximations successfully 
follow the contortions of Vs, but as t increases V\ loses accuracy. The initial dot labels V^O), the second dot 
corresponds to V$(2), and the arrows label Vs, V\ and V 2 at t — 22.5. The second order approximation V 2 holds on 
longer, but eventually succumbs. By then, Vs is far from the constant Lie quadratic Vq. □ 



Theorems |2j [3j give explicit formulae in terms of elementary functions for approximations V n to nearly constant 
Lie quadratics Vs in E 3 . Order n approximations x to Riemannian cubics xs in SO(3) with Lie quadratic Vs can 
be found by solving the linear system of ODEs 

z«(t) = x(t)V n (t) (16) 



Higher order Taylor polynomials are almost as uncompetitive. 



7 




8 




Figure 2. V s (t) (blue), Vi(t) (green), %(t) (red) for t e [0,25] 
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where, on the right hand side, V n (t) is identified with an element of so(3) and, considered as a 3 x 3 matrix, 
is premultiplied by the unknown matrix x(t). Because (16) has nonconstant coefficients, x cannot be written 



down directly and would be found by a numerical integrator such as NDSolve. Considering that Riemannian 
cubics in 5*0(3) are solutions of the 36 dimensional nonlinear system |lj), the 9 dimensional linear system (16) 
appears comparatively benign. But solving ( 16 ) is nontrivial, and this step can be avoided: an explicit first order 
approximation x to xg is given as follows. 

6 Approximating Nearly Geodesic Cubics in 5*0(3) 

Given a generic Lie quadratic in a semisimple Lie algebra, there is an integrability algorithm |13j that takes a 
single quadrature, and gives an explicit formula for an associated Riemannian cubic. The algorithm takes a simple 
form |12j for Riemannian cubics in 5*0(3). So we might try to approximate a nearly geodesic Riemannian cubic 
xg : [t Q , t-y] — > 50(3) to first order by substituting V% for Vg in the integrability algorithm. This has the unexpected 
benefit of giving a first order approximation for the quadrature in terms of elementary functions. On the other 
hand, V% is needed to approximate the other terms in the integrability algorithm to first order. This gives, by 
somewhat indirect means, explicit first order approximations x to nearly geodesic Riemannian cubics x in 50(3), 
without even the need for a single quadrature. To review the integrability algorithm for Lie quadratics we make 
two definitions. 

Definition 2 For r € K, let R r £ 50(3) be clockwise rotation by r in E 2 , namely 



cos r sin r 
— sin r cos r 
1 



□ 



Definition 3 For linearly independent X\,X 2 S E 3 , set 

\\Xl \\X 2 — (Xi ,X 2 )Xi 



S(X U X 2 ) 



X ± xX 2 



^/\\x 1 p\\x 2 p-(x 1 ,x 2 y 2 ^/\\x 1 p\\x 2 p~(x 1 ,x 2 y 



X! 

IIXill 



i T 



e 50(3) 



where T means transpose. Then, for any curve W : [io^i] E 3 with W{t), W^\t) everywhere linearly indepen- 
dent, define T{W) : [i ,£i] 50(3) by 



T{W)(t) := S(W(t),W {1) (t)). 



□ 



Now we return to Lie quadratic^W . As is easily proved [TU] , they have constant acceleration: c :— {V^ (t) , (<)) 
is constant as t varies. When 7™ vanishes identically in a nonempty open interval, the corresponding Riemannian 
cubics are cubically reparameterised geodesies. Assume V^(t) ^ for any t e [tcb^i]; where V : [toj*i] ~^ so(3) = 
E 3 , and E 3 is identified with so(3) in the standard way, by the adjoint reprsentation. 

Lemma 1 Let x : [to,ii] —> 50(3) be a Riemannian cubic whose Lie quadratic is V. For t G [^Oi^]; s ^t 



y(t) := R m oT{V^){t) where 0(t) := c 1 ^ 
and o stands for matrix multiplication. Then x(t) = x(to)y(to) T y(t) . 



•-(c,v (2) W) 

\\V( 3 )(s)\\ 2 



Proof: Since left Lie reductions are invariant with respect to left multiplication, we can suppose without loss 
of generality that x(t ) is the identity. Setting W 3 (t) := c^^V^ (t), Wi{t) := V^^/WV^^W, and W 2 (t) := 
Ws(t) x Wi(t), we find that W\,W 2 , W3 meet the requirements of Theorem 5 of [12] (in [12] a precise choice of W\ 



3 Here the symbol V is used to denote a Lie quadratic, rather than a variation of Lie quadratics. 
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is not made, allowing other curves in S 2 orthogonal to W3). So there is a Riemannian cubic y of the form above, 
with Lie quadratic V, provided 

^(t) = (W 1 (t),W^ ) (t) + nd- 1 (V(t))xW 2 (t)) = -(w[ 1) (t) + ^ 1 (V(t))xW 1 (t),W 2 (t)). 

With the present choice of Wi, = . r ,„, lo (V (3) , [y (2) ,F (4) ] + [V, \V {2 \ V {3) }}) = 

c I \\v v d ) M 



(yO), [yC2), [y(3), v] + pA 2 ), + [V, [T/< 2 \ V^}}) (V^, [V™, [V^M 



C l/2||^(3)||2 c l/2|jy(3)||2 

u, (V^\V^) = r U2 {V {2 \[V {1 \V]) = m(H^>) 

' ||v( 3 )|| 2 ' ll^ (3) ll 2 ' ll^ (3) ll 2 

□ 

So Riemannian cubics in 50(3) can be found from Lie quadratics by a single quadrature. First order approximations 
to nearly geodesic Riemannian cubics can be written explicitly in terms of elementary functions, without the need for 
quadrature. More precisely, let x$ : [to, t±] — > 50(3) be a Riemannian cubic whose Lie quadratic Vs : [iojii] — > so(3) 
is nearly constant. For n = 1, 2 let y„ : [t , t±] — > so(3) be the order n approximations to Vs given in q5l Assuming 
B^O, write B = /3(cos7, sin 7) where 7 € [0,2vr). Set p := -2c 2 /(d 2 /3). 

Theorem 4 For i e io 5 *i] define x(t) :— xs(to)y(t ) T y(t) where y(t) := R^u\ T(V 2 )(t) and 

•= ^( p 2 , i)i/2(( f 1 aii(cos(7-^-^o))-cos7)+Qi 2 (sin(7-rf(f-fo))-sin7) 

T/ien /or j > 0, and all t G [to>*i]> x s (*) = zc^(t) + 0(S 2 ), where the asymptotic constants depend on 
out ai2) /3j c 2,d, to>*i ondj. 

a A (2) * (3) 

Proof: Let W^t) and Wi(t) be the unit vectors in the directions of V 2 (*) an d ^ (t) respectively. Taking 

V = V s in Lemma [l] it suffices to show that cj)(t) = 4>(t) + 0{5 2 ) and that W t (t) = W t (t) + 0(6 2 ) for % = 3, 1. 

By Theorem g v/ 2) = ^ (2) + 0(S 2 ) = v[ 2) (l + O(S)). So c = c(l + 0(<5)) and C = 6(1 + 0(5)) where 
c : = \\v} 2) \\ 2 = <S 2 (4c 2 + d 4 /3 2 ), C := <5 ad(2c 2 , -da 12 , da xl ) = ^ (2) - fc (1) ,Vi] + 0(£ 2 ). Similarly y/ 3) = 
Vi (3) (l + 0(5)) =>> ||V5 (3) || 2 = rf 6 /3 2 + 0((5). So 

_ 1/2 /' c-(c,^ (2) ( s )) _ /2 c-(c,v; (2) ( s )) 



Because the denominator is constant, the integral on the right can be computed precisely, giving 

<j)(t) = 5(ic 2 I rf 4 /3 2 ) 1/2 ( t ~ t ° I a ii( cos (T , -o'( f -*o)) -cos 7) + Q 12 (sin(7-o;ft-t )) -sin 7) ^ | 2 
2 d 2 d 4 j3 

= 4>(t) + 0(8 2 ). 

By Theorem^ V s {2} = V 2 {2) + 0(S 3 ) = V 2 {2) (1 + 0(5 2 )). So W 3 = W 3 + 0(S 2 ). Similarly W x = W x + 0(S 2 ). □ 

Whereas in [11] the analysis of Riemannian cubics in 50(3) is complicated, in the present paper x is algebraic in 
terms of polynomials and trigonometric functions, as we see by taking n = 2 in equations (|9]): 

/ 2 (2) = 2d[A ,i(e-l)B] + [A l ,2(e-l + d(t-t Q )ie)B] 

v { 2 2) = -^(i-l)A Q + ^(d(t-t )l-Xe-l))A 1 +2(2c 2 (t-to)+d 2 I(q))t£B 
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where everything is constant, except t, 
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cos{d(t - t )) sm(d{t - t )) 
-sin(d(*-i )) cos(d(*-* )) 



and I(q)(t) 



C2{t-t Q ) 3 , ci(t-to) 5 



+ c (t - t Q ). 



Then V 2 {2) = S(2c 2 F - d 2 £B) + f (/ 2 (2) F + v [ 2 ] ) is affinc in e(t)B, / 2 (2) (t) and i^W- So in TheoremEJ T(F 2 (2) ) 



is 5(X 1 ,X 2 ) where 



X 1 (t) := 2c 2 Fo~d 2 £(t)B+-(f?\t)Fo + vi 2 \t)) and 



d 3 ie(t)B+-(A 3) (t) + v { 2 3 \tj). 



Example 10 For some small 5 € R, Zet V5 : [0, 10] — > so(3) 6e i/ie Lie quadratic satisfying 

ad(Vi(0)) = (1,0,0) +5(0, 1,0), ad(yj (1) (0)) = 5(0,0,l)/2, ad(Vi(2)(0)) = 5(l,l,l)/4 (17) 
and Zei Xg : [0,8] — > 50(3) be the corresponding Riemannian cubic for which Xg(0) = 1. 

Given any particular value of 5, say S = 0.05, Mathematica's NDSolve can be used to numerically solve the 
quadratic differential equation |7p for Vs : [0,8] — > so(3). Then x$ is found by numerically solving the linear 
differential equation 

xf\t) = x s (t)V 5 (t) 

where the coefficients on the right are entries of the matrix V$(t). The second rows of xs(t) are shown as the blue 
curve in Figure^ with blue labelled points corresponding to t = 0, 1, 2, . . . 10. A geodesic represented in this way 
would appear as a circular arc. 

To compute the approximate solution x : [0, 10] — > SO (3) from Theorem^ take to = 0, t\ = 10 and D = (1,0,0). 
For x with x(0) = 1 to satisfy the other initial conditions (i7[ ) with 0{8 2 ) errors, it suffices to have 

Vi(0)) = (1,5,0), Vf>(0) = (0,0,5/2), Vi(2)(0) = 5(l,l,l)/4 

where V\ is considered as a curve in E 3 . For this, choose parameters Co, c\, c 2 , aoi> Q-02, an, ai2, ft, 7 so that 

(c ,a i +/3cos7,a 02 +(3sm-f) = (0,1,0) 
(ci,an + /3sin7, a i2 - /?cos7) = (0,0, l)/2 
(c 2 ,-/3 cos 7,-/3 sin 7) = (l,l,l)/4 

Taking (i = \2/A we have Co = C\ — 0, c 2 = 1/8, 7 = 57r/4, aoi = 5/4, a 02 = 1/4, an = a 12 = 1/4. The 
second rows of x(t) are shown as the red curve in Figure [3j with red points corresponding to t = 3, 4, 5, . . . 10. For 
t = 0,1,2, x(t) and x$(t) are hard to distinguish. For t = 3,4,5 the approximation is good enough for the blue 
numbers to label both. 

It is hard to see much difference between x(t) and Xg(t) until around t = 4, and long before then the Riemannian 
cubic x$ is obviously nongeodesic. After that, differences between x(t) and xg(t) become noticeable (t — 5), and 
then large (t > 6) where x$ is very far from geodesic. □ 

7 Conclusion 

We derive first order approximations x to nearly geodesic Riemannian cubics xg in 50(3) defined over a given 
interval [to, ti]. The approximations capture a good deal of the geometry of xg, including behaviour of the associated 
derivative, corresponding to the body angular velocity for a rigid body. Equivalently, the approximations apply 
to restrictions of arbitrary Lie quadratics to sufficiently small intervals. This fills a gap in the literature, by 
complementing studies on long term asymptotics of Riemannian cubics. 
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Figure 3. Approximation (second rows) to a Nearly Geodesic Riemannian Cubic in 50(3) 
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